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ABSTRACT 


Certain experiments contemplated for space platforms must be isolated from 
the accelerations of the platform. In this paper an optimal active control is 
developed for microgravity vibration isolation, using constant state feedback gains 
(identical to those obtained from the Linear Quadratic Regulator [LQRj appnach) 
along with constant feedforward (preview) gains. 

The quadratic cost, function for this control algorithm effectively weights 
external accelerations of the platform disturbances by a factor proportional to 
(1/a.')*. Low frequency accelerations (less than 50 Hz) an' attenuated by gr< atei 
than two orders of magnitude. The control relies on the absolute position and 
velocity feedback of the ex; eriment and the absolute position and velocity 
fml— forward of the platform, and generally derives the stability robustness 
characteristics guarant<vd by the LQR approach to optimality. 

Tht' method as derived is extendable to the case 1 in which only the relative 
positions and velocities and the absolute accelerations of the experiment and space 
platform are available. 



1. INTRODUCTION 

A space platform experiences local, low frequency accelerations <0.01-30 Hz) 
due to equipment motions and vibrations, and to crew activity [l]. Certain 
experiments, such as the growth of isotropic crystals, require an environment in 
which the accelerations amount to only a few micro-g's [2]. Such an environment is 

not presently available on manned space platforms. 

Since the experiment, and space platform centers of gravit y do not coincide, a 
means is needed to prevent the experiment from drifting into its own orbital motion 
and into the space plat form wall. Additionally, some experiments require umbilicals 
to provide power, experiment, control, coolant flow, communit ations linkage, ot 
other services. Unfortunately, such measures also mean that unwanted platlorm 
accelerations will In* transmitted to the experiments. This necessitates experiment 
isolation. Passive isolators, however, cannot, compensate for umbilical stillness, nor 
can they achieve low enough corner frequencies even if umbilicals are absent. 
Active isolation is therefore essential. 

The problem, then, is to design an active* isolation system to minimize these 
undesired acceleration transmissions, while achieving adequate stability maigins and 
system robustness. Spatial and control energy limitations must also be 
accommodated. 
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2. MATHEMATICAL MODEL 

Tho general problem has three translational and three rotational decrees of 
freedom. For simplicity, however, this analysis will consider only the 
one— dimensional problem. The general problem could be treated in an analogous 
manner. Let the experiment be modeled as a mass m. with position x(t). Assume 
that the space station has position d(t), and that umbilicals with stiffness k and 
damping c connect the experiment and space station. Suppose further that a 
magnetic actuator applies a control force proportional to the applied current i(tb 
with proportionality constant o. Such a model is shown in 1 igmc L 
The system equation of motion is 

nix + c( x — <1 ) + k(x— d) + oi = 0 


Division by m and rearrattgement yields 


x 



(N-<I) — ^ (x— t'O - 
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In state space notation this becomes 


x = A x + b u + f 


where 
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f = 


— (i + £ (1 
in m 


The objective is to minimize the acceleration x(t). 
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OPTIMAL CONTROL PROBLEM 

The optimal control problem is that of determining the onuol current 
u(t) = i which minimizes a suitable performance index 


J = J(x, u.t) l u 

for the system described by Eqii . (3) subject to the state variable* conditions 

x(0) -- x Oat 

1 i m x(t ) --- 0 b>bi 

t — X 

Another reasonabh' assumption is that f(t) is l>ou tided. and it will be hand 
mathematically advantageous (and only minimally restrictive*) to assume that [! t ) t- 
also a dwindling lunci ion 

I i m f(t) = 0 (:>r ‘ 

t — x 


A quadratic performance index 

•1 = -yf [x 1 VV x + w .11“] dt C (l 

“0 

has been chosen, as one that lends itself well to the variational approach to optimal 
controls, since an analytical solution is desired. The upper limit of the detiuite 
integral has been selected so its to yield a time-invariant controller Here Wj is a 
square 2x2 constant weighting matrix while vv.j is a weighting constant. 



Although, W, could be a full 2x2 matrix, for this problem a diagonal lorm 
has boon employed lor the sake of simplicity. 


Wf = 


vv 


la 
0 w 


0 

lb 


(T 


The performance index consequently reduces to 


.) 



w 1 b x 2 



so that each state is weighted independently. 

If sinusoidal motion of the experiment is considered, so that 


x(t) = B sin 0 , 1 . 


and x(t) 


J- X (t), the cost function can be expressed in t« 


ms of the acceleration 


and cont rol as 


X vv 

l(- 

0 


i W . i ■) •) *) 

U + -4]) B"x“ + w.^if | dt 


It is apparent that this performance index couv< 
frequencies much more tin. at higher frequencies. 


*i den i lv weights acceletio ion • ,0 a>. 
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I. SOLCTION 


Finding the optimal control to minimize Eqn. (4) is a variational problem of 
Lagrange, for which the initial steps of the solution are well-known fe.g.. Elbert 
[4]). The variational approach is outlined Itelow, following which the complicat ions 
added by the non homogeneous term f(t) will be addressed. Current optimal 
controls texts either assume that f ( t ) = 0 (e.g., [4], p. 2f>2) or require that it have a 
restricted range space (e.g., [6], p. 238). The solution that follows provides an 
analytical optimal without imposing such restrictions. 

The argument of the cost function J from Eqn. (4) is augmented by the 
Lagrange multiplier A times the system equation of motion Eqn. (3) where 


The result .1 can be expressed as 

•X’ 


II dt 


(!!■ 


0 


where the Hamiltonian 11 is 


11 = ^ (x 1 W j x + w.jiC) + A^ (x — Ax — bu 


It is desired to obtain an optimal solution u = u which minimizes .J. 
The first, variation of J(x, u, x) is 


wind, is sot, equal to w*n> to minimize J. 


llow<'V(T, integrating by part s, 


rOO 
• 0 


(^1 fik) dt, = 
i)k 



/>X (it 


so that the ab >ve expression for M becomes 


,-oc 

6J = 

j 0 



(13) 


Bot h -Sx and bn are arbitrary variations, so M - 0 only 


0 \\ _ \ 1 
lft~ ~ 


_ 

7ft- 


1 Ibi 


The conditions given by Eqn. (5) sli11 a PP 1 . v - 
Solving Eqns. (Ha) and (Mb) yields 


A = VV j x - A A 


(15a) 


u 


w 


I ) 1 A 


(15b 
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* 

Temporarily eliminating u produces the result 


fxl - fxl 


M 

7 = A 

• + • 


UJ IaJ 

i 

O ! 


(10) 


whoro 
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If Kqn. ( 1 (i) is now solved lor A in lor ms ot x ami [. If j u . (15b) will i hen I m mi.di an 

f 

expression for the optimal control w . 

As noted before, optimal control texts generally treat the homogeneous 
problem (where f(t) - 0), but. they do not provide an analytical solution to the 
non homogeneous system described by (5) and (lb). Salukvadze has treat d the 
non homogeneous problem [•1,5], but his dilficult. treatment seems largely ! 1 have 
remained either uncoinprehended or under— appreciated. 1 his me - hod is especially 
well-suited to low-frequency disturbance rejection, and lias been applied below <o 
the present problem. 

'The homogeneous solution to Eqn. (15), where t = 0, is 
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The four eigenvalues of A may be found to be, in ascending order of real parts. 




lb = ~ 


[ -i\ + (/V 2 - 4 ^2 )1/2 1 
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[ ~l\ - U\- - 4/^) 
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fl :] ~ /4 l 

/* j ~ — / 4 2 


(18a) 


( 18b) 


(18d) 


where and ,j.> are defined as follows: 


and 
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(19a) 


(19b) 


The eigenvectors of A corresponding to the respective eigenvalues /i k may be chosen 
to be 

l 

"k 

7. t 7[ 2 7,(7 2 + /' k ) . 

— + — — + 

/‘k Vk 

?1 + ^ 7 o + I 1 k )/'k 


Hk = 


(20a) 


where 7^ , 7,-,, 7^, and 7^ are defined below: 
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C20b) 
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(20c) 


( 20 ( 1 ) 
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7. = w, 

1 la 


(200) 


I'sing Eqns. ( 18 ) through ( 20 ) with ( 17 ) the solution to the homogeneous system i.- 
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K J 


with m = j ''l J- , k = 1 , 1 ami where ( j are arbitrary constants. 


Application of the variation of parameters method with terminal conditions 
(Eqns. 5b. c) leads to the general solution of the non— homogeneous system, with two 
constants of integration yet undetermined. 

If the two constants of integration are eliminated by solving for A in terms of 
x and f, the general solutions for Aj and \. y become: 


Aj — + C) x 2 + £-jO 


-/‘it 
■ + 

c ( > 

s 4 ( 

(22a) 
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(22b) 
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in which the ^'s are functions of the eigenvalues and eigenvectors ol A, and of the 
disturbance f(t). 

Thp Solution Form 

Using the fact that 

u * (t) = i_ A T b [cf. Eqn. (15b)] (23) 


the optimal control is found to be 


u (t) = T/jXj + 

— //,t fit t l l 2 l 

+ /y.jC / e t.j(i )<lt + V,f / 0 ( 

(24a) 

where 

'll = =7T 

(24b) 


1-> = ^ b + '‘I + ,l i ] 

(21c) 



(2 Id) 


~ 77 ^/t j ^ /( 2 + m 1 2 m ^ 

(24e) 

(It should be noted that the feedback gains ii { and are those wIikh vu,llltl 

result from applying standard LQR theory to the homoge.teous system equation 
x = Ax + bu). In Eqns. (24) /x., are the eigenvalues of A with negative real 


parts, [see Eqns. (18a,b)| and 
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By repeated application of the method of integration by parts, the control may be 
re-expressed in terms of an infinite sum: 


ti 0) = Vi + 'h*2 + ,; 3 


oo 
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r = 0 


(-1)’ f., (r) (t) 
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r + 1 


+ 7 4 


foe (-1)' r.,( r >(t)l 
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^ n r+1 I 

r = () //. 2 J 


(25) 


Rewriting f., in terms of d and d. t he control function becomes 


II rn = //jX(t) + /A_,x(t.) + + ~)} d(t ) 


r— 1 
+ 'J 
i-l 


>i-l c , , + 


a k , 


■It ' s; (— f + ~f) + i‘ m t?t + ~ rlr’ 

/< i N 


m 1 l + l 

th lb 


d (i) (t. 


+ |( — i ) n * J_ ( — ^ — I — )j cl^ 11 \t) + higher order terms 

ft i 11 lb" 

1 “ ( 26 ) 


This may be written in a more appealing form as 


u (t.) = c () x(t) + < v x(t) + c (|0 d(t ) + c (il 


d(t) + higher order terms 


in which the constant coefficients e , c v , c ( jy, and c. ( jj may be defined from Kepis 
(24) and (26). Clearly, if the infinite sums converge rapidly enough, the optimal 
control can be approximated by 



u*(t) = x(t) + c v x(t) + c d0 (l(t) + c (1 , (1C ) 


(28) 


For very low frequency disturbances the higher order terms in Kqn. (26) are 
negligibly small, and the control (Eqn. (2«)) closely approximates the optimal. If. 
in fact, the second- and higher-order derivatives of d(t) are identically zero, the 
approximation is exact. It can be shown that for the critically damped closed loop 
system the eigenvalues are real and equal, and that the convergence is more rapid 
than for the overdamped system. Further, as the closed-loop system eigenvalues 
become more negative the convergence speed goes up as well. 
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5. CONTROL EVALUATION 
Physical Realizability of the Control 


The control, Eqn. (25), is physically realizable, if the states and sufficien' 
derivatives of d(t) are accessible (or estimable by an observer), and if the higher 
order terms are negligible. It is not necessary that the eigenvalues be real, although 
the proof of this requires a more general linear-algebra or state-transition-matrix 
approach. 

If values an* assigned to the system parameters, associated controller gains 
can be evaluated. Suppose that m = 100 Ibm, k = 0.3 Ibf/ft, c = 0 Ibt— sec/lt. and 
o = L0 lbf/Amp. With w.^ arbitrarily set at 1 and Wj () vari »d. associated integer 
values of Wj can be found Itelow which the eigenvalues /tj mid //., will always be 
real. Such values are tabulated in Table 1. Stated otherwise, the tabulated values 
of the weights Wj and w (l) are those integer values (for the mke of simpliciiv) for 
which the closed loop system is closest to being critically d.unped without being 
underdamped. Corresponding controller feedback and feed-forward gains (for the 
first five derivatives) tire also included. 

The states x(t) and x(t) and the derivatives d (n '(t). d’ "ltd and d"'(t) are 
clearly available for an earth-based system. However, in space, me only absolute 

measurements which can be directly available are x(t) and d(t). Irom which xtt). 
c'l(t) and x(t), d(t) are obtainable only by successive integration(s). Rearrangement 
of (2S) into 


.T(t) = (c p + I ||(J )S(I) + (C v + >, n )N(t) -C ( |0 |x(n -.1(1 I , -ri(t)] 

( 20 ! 

or 

u*(t) = (c + c (10 )d(t) + (c v + c ( 1 J )d(t) + C p [x(t) -d(t)] + c v [x(t) -3(1.)] 

( 30 ) 
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obviates the need for one accelerometer, but one accelerometer plus nvo integrations 
remain necessary for either the platform or the experiment. Since |x(t>— <lft >1 (or 
one of its integrals) has not been weighted in the performance index J. experiment 
drift will be a problem that must be corrected either by another control loop or by a 
change of system states. The latter could be accomplished by incorporating an 
accelerometer attached to the experiment into the state equation. Alternatively. 

one could append legrator to the plant, include the current i(t) as a third slat,-. 

and optimise the control di/dt. But lor the sake ol simplicity B f i ’si 


former has been assumed (without developn I in this |>aper. 

The higher order terms of tire control [Eqns. (Ml and 2(1)1 can is- neglected, 
for low frequencies, if the eigenvalues p [ and //., are of surfretent modulus. 1 " 

eigenvalues. in turn, rontrol of the der igner. determined l.y Iris choice 

of weights w u , w M> , and w.,. it is apparent from '.., 1 ,. (in) Ural u 

reduces to iwo alternating power series, l or a sinusoidal disturbance of frequency - 
tire series form of the control converges tor I a//c | <1 (i - I — I- I' 1 an 

1 hat each alternating power series converges like ^ ^ (-1) I (l ) 

frequency rlistnrhances (i.e., small relative u, system dosed loop eigenvalues) a 
™,„rol formed by series truncation very closely approximates lire optimal- 

For example, suppose that the ruralized frequencies U7/<;! for n sinusoidal 

disturbance are less Ilian 1/5, and Ural only the feedforward control terms <,, () dir! 
ami c ()1 d(t) arc inclmhxl with the feedback terms. Even so. the Mforuard porlion 
of the' truncated control, at any time t. will be a current that is still within 1 • [>•■■„ 
(1/5)" of the feedforward isrrtiou rrl the actual optimal. II llu 
frequencies are below 1/10, this approximation error will Ire less than U. table 1 
shows that the gains c (|j of higher order derivatives d" l 0) (see Eqn. (hi) lor 
algebraic representations] are, in tact, quite small. 


q r; h ^ pc 
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In some circii instances there may be design constraints which prevent the 
designer from selecting weights that will lead to sufficiently rapid convergence. 
However, since convergence occurs rapidlv even for eigenvalues of relatively small 
modulus ( | 1 < 1/:?), in a great many cases the designer will have much 

latitude in his choice of weights. For "low" frequency disturbances, in these cases, a 
control which includes only one or two feedforward terms will be "close" to the 
optimal. These frequencies will be well— attenuated. 

Higher frequency disturbances will also be well— attenuated, provided the 
input— to-output transfer function(s) are at least strictly proper in the I.apiace 
Transform variable s. This will not be the case for the present problem if more than 
three feedforward gains (c^, c^, c^.,) are included in the coin ml. Practically, this 
means that only proportional and first- derivative feedforward [Eqn. f 25 j with r = 
0.1 or Eqn. (26) with n = 2j should be added to the feedback control terms. As will 
be seen shortly, however, adding even the proportional feedlorward term(s) can 
dramatically improve the disturbance rejection over that afforded by LQR feedback 
alone. 


Transfer Function and Block Diagram 

Neglecting the higher order terms, the transfer function Ix'tween input and 
output accelerations or displacements is 
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and a block diagram of the controlled system can be drawn as in Figure 2. 



Control Stability. Stability Robus tness, and General Robustness 

Since the control feedback gains are the same as those obtained by solution of 
the standard Linear Quadratic Regulator (LQR) problem, the closed loop system is 
stable and enjoys the stability robustness characteristics guaranteed by the (LQR) 
approach to optimality, viz., a minimum of 60° phase margin, infinite positive gam 


margin, and 6 dB negative gain margin [6], Additionally, numerical checks indicate 
that it enjoys substantial insensitivity, or general robustness to uncertainties in k, c. 
and m, as indicated by Table 2 and Figures 3 through 10. By comparing the Bode 
plots of Figures 3, 5, 7, and 9 (corresponding to controls using both LQR F/B and 
proportional F/F) with those of Figures 4, 6, 8, and 10. respectively (corresponding 
to controls using LQR F/B only), one can see that adding feed-forward 
substantially improves disturbance rejection at low frequencies. For example a 
comparison of Figures 3 with Figure 1 indicates that the* optimal control method 
described above can lead to acceleration reductions of greater than tour orders oi 
magnitude for all frequencies. This reduction is more than two orders of magnitude 
below that afforded by LQR feedback alone at the lower frequencies, i.e., those most 

heavily weighted in the performance index. 

The order of the reduction is eventually limited by control cost, ol course, 
probably in terms either of actuator— related limitations (such as hint icmoval 
force-generation requirements) or of power limitations (especially in a space-station 
environment). The control also leads to displacement reductions of the same 
magnitude, limited in this case by actuator-stroke or spatial limitations. Providing 
a unit transmissibility for very low frequencies and weighting (x-d) and/or /(x-d) 


in the performance index .1 
limitations. 


would be steps toward addressing these lattei 
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Computational Aspects 

A significant amount of algebra was required to solve the two— state problem 
of this paper, and the labor involved increases dramatically with each additional 
state. However, such symbolic manipulators as MACSYMA may be used to ease 
the workload if a symbolic solution is desired. Further, well-known numerical 
methods exist (i.e., Potter's method [7] or Laub's method [8]) for solving the 
solution to the homogeneous system. These can readily provide the feedback gains 
in numerical form, even for problems with many states. It might be anticipated, 
then, that a numerical method also exists for finding the desire 1 feed— forward gains. 
Such is the case, as will l»e shown in a later paper. 



6. CONCLUSIONS 

This paper has applied an existing method for obtaining an optimal control 
to the microgravity platform isolation problem, for which the disturbances to be 
rejected are low-frequency accelerations. The system was assumed to be 

representable in the form x = Ax + bu + f, with quadratic cost function 

j _ 1 °°(x T W.x + w,,u 2 )dt and diagonal weighting matrix W’y l'he. resultant 

-Jo 1, 

control law was found to be simple, stable, robust, and physically realizable. 
Further it was shown to have excellent acceleration- and displacement-attenuation 
characteristics, and to be frequency— weighted toward the low end of the acceleration 
spectrum. 

'l'he method is extendable to the case for which only relative positions and 


velocities, and absolute accelerations, are available; and can be applied so as to 

weight relative displacements in t he poi'lormanee index. 

The approach as presented is algebraically intensive, but symbolic 
manipulators can be used to ease the algebraic .abors. Further, since the method 
produces feedback gains identical to those obtained by the LQR approach to 
optimality, numerical computation of those gains is easily accomplished, even lot 
large systems. The feed-forward gains can be found numerically with comparable 


ease. 
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Table 1. Optimal F/F and F/B Gains for Selected State 
Variable and Control Weightings. 

System Parameters: 

m = 100 Ibm k = 0.5 lbf/ft 

c = 0.000622 lbf-see/ft (< = 0.17.) 
n = 10 lbf/amp 
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Table 2. Closed loop transfer functions for system with design 

parameter values of k = 0.3. c = 0.000622, and m = 100; but 
with actual parameter values as shown. Gl, G3. 65. and Gi 
include both LQR F/B and proportional F/F; 62, 64, 66. and 63 
include LQR F/B alone. Weighting parameters used were 
wj = 253, w lb = 10, w 3 = 1 (see Table l). 


System Parameters 


/ lbf \ 
k ( fr ] 

, ,lbf-sec, 

c * Ft > 

m (lbm) 

0.3 

0.000622 

(C=o.l7.) 

100 

0.3 

0.000622 

100 

0.45 

0.000622 

100 

0.15 

0.000622 

100 

0.3 

0.00622 

100 

0.3 

0.00622 

100 

0. 15 

0.00622 

90 

0.15 

0.00622 

90 


Closed Loop Transfer Function 

* 2 X(s) 

s-D(s) 




6 1 ( s ) 

0.0000622s + 0.0001 ! 

- 3 . 

0 . 31056s"f4 . 4675s + 16. 0624 

6 2 ( s ) 

0.0000622s + 0.0300 
0 . 31056s 2 +4 . 4675.S+ 1 0 . 062 1 

G3(s) 

0.0000622s + 0.0151 
0.31056s* +4 ,4675s+ 16.077 1 

G4(s) 

0.0000622s + 0.0450 
0. :J105f)S-+4. 46 |5S+ 16. Oi - l 

65 (s) 

0.000622s + 0.0001 
0 . 31056 s 2 +4 . 4680s+ 16 . 0624 j 

C6(s) 

0.000622s + 0.0300 
0 . 31056s 2 +4 . 4630s+ 16 . 062 1 j 


, 0.000622s + O.OIol 

6 1 ( s ) = r; 

0. 27950s" +4 . 4680s + 16.01 . 1 


, 0.000622s + 0.01)0 

6<3(S) = 7 — ' 

0. 27950s"- l.lOSOsr 10.0. .4 
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DISTURBANCE LEVELS 



Periodic Accelerations 


Re' ' 


orijlii 


2E-2 


2E-3 


enc " 


5 to 20 


Source 


Thruster Fire 
(orbital) 

Crew Motion 


2E-4 


Ku Band Antenna 


Non-Periodic Accelerations 


Relative Gravi 


IE-4 


Source 


Thruster Fire 
i (Attitudina D 


IE-4 


1 


Crew Push-Off 
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General Probl em : 

In the upcoming space 
station, planned by NASA 

for completion in trbe |^50s, 
minimize the low freauenc 


accelerations transmitted 

from the space station to an 
experimental platform contained 
on (inside) the space station. 

“minimize"; reduce to ~ 10 % 

• L I 

it possible 

“low frequency": o.ooi to 20 Hz 
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x x x_x x x 


x x 


x x x_x x x 


x x x x 


x x 


x x H'K X x x X m X X 

x xjj X X x X X X 

The Lorentz Equation: F = ilXB 
F = Force 1 = length of wire i = current 

X Represents the tail feathers of a magnetic field B vector into the page 
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Objective: 


Find tbe "be st* % i (*), 
to minimize x (+) . 
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State Equations: 


o I 
m m 


+ -j* t+ k d+ £d 

I m I ( m m 


O ft 


b u 


A = [-I -il 

L m m J 

>■[•;] 

‘ ■ ■ 1 4 1 ■ ' U4 


456 



University of Virginia 


ROMAC 


Problem Statement : 


Determine the control u(t) 

which minimizes the performance 
index £ 

f 1 ( x T W, x + u T W u)dt 


for 


$u 


the system 

x * A * ♦ 3 u 


Jjject to the conditions 

J x (o) - X o 

I I im X (t ) — O 


I i m f (t) * O 


t-+ 
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Solution Method ( D i ff e rent- i a I 

Equations Approach) : 

•• Augment the performance index J 
with the state equations using 
Lagrange multipliers. 

2. Take the I s * variation S J of -the 
aujwen + ed performance index J and 
Set it e^ual to zero: 

s J® jf [ir ^5 + "It Su + Sxj<it - o 

where H= f (x T W, x + u T W, u) + * T (s- A X-B u - f) 

3. Integrate the third term of ■the integrand 
by partSj combine terms^ and set coefficients 
of the arbitrary variations Sx and Su 
ec^ual to zero. 

Result; A - W| X - A A 

a - w/bV 

t. Substitute for U in the state equations, 

to yield 
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tk. solution ({.«, the .pti—i ts = y*) 

is n6 w u* = W,-'B T X where A it *•««* b y 
Solving the *ys+em 

( i\ [ft w 3 -V](<K( f , 

W “|.W, -A T JhJ \°-l 

subject to 

X Co) a X* 

lim xW = 2 
*-* *° 

lim £ £t) a 2 

+ -♦•0 

s Find the solution of the homogeneous system 

,.*(*).[* bw;'b t fa] 

’ Ui k -* T JUi 

t . Use the variation of parameters method 
to find the general solution of the 
nonhomogeneou. system, i.e., of 

(i\ A BW 3 "B t ] fS ] + (f 

ark, -a t JUI ) 

7 Apply the terminal conditions on 
' conclude that n of the Zn arbitrary 
constant, in the general solut.on are 
* n Lial to zero. 
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9. Sftlvc for (t) if* 'tcrwi of X (t) m a manner 
S t4 c h that *tKc rcmdin i n orb i ^ retry 
constant# are eliminated. 

Use the e^u at i on 

a*- W,-'b t > 

to ■Find u* in terms of X. 

Result; 


u*w « (Wj * 1 B T x tl x;')x 

♦« » T C) 




where 


0 - A 


is the Tordan 


Canonical Form of the HamiltosioJ 

matrix Ta 

J ) 

w. -a t 

U 1 J 

where A contains only the negative 
eigenvalues of the Hamiltonian 
matrix^ corresponding to the 
eigenvalues of the closed -loop 
System (assuming ^ contro liable^ 
where X s X# X fl the eigenvector 

,X t | X 2lj 

matrix which leads to t ne 
above J. C. F. ^ 
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c-l) 

xr 

II 

n 

(-•) 

v c-rt 

11 

An 

ion 

of 


-j 

integral requires constants 
of integrati on that a re all 
identically zero. 

to. Integrating repeatedly by parts and using 
the •fact* that 

x::yx;r= x.ax; = a= a- bw,-b t p 

where Pd'-X^X,*' [? i* the solution to the 

well- known Algebraic 
Riccati Equation] j 

develop equivalent forms for u *60 : 

u» ( t) B -vo t p* + w;'B T x;rv A ^x;>)dt 

, - W ," 1 B t P x - W,-B T X^*' r I 


= - w:' B T P X - Wj' B T z(-A" T ) r "" p f <rt 

3 J roO 

= - w,"' B T P x - %' B T x; T r i(-A'T-x; f - 

ote-. A state transition matrix approach yields 

a *M . -w," b t p S («‘X>dr 


46 




University of Virginia 


ROMAC 


S o I u f i o n » 

u‘(+) = -W 3 -'B t P x + Wj" B T X‘;“"e A+ f e A * dt 

= - w:' b t p x - w:‘ b t x‘;°V a * f t AT x‘; rt f m ** 

3 - 3 21 ^ II - 


= -w ; 1 b t p x - w ; 1 b t x;r 1 (-A-r r, x:;''f 

3 - 3 11 ris0 11 


#~i »\T »/ (- 0-1 


^,vr + l v C -0 _<r) 


= -w 3 *'b t p x - w:‘ b t z (-A" T ) r+, Pf <r * 

r» o 

= - w;‘ b t p x - w 3 “ b t x; T £ xj f <r) 


Dropping higher order "terms (for r > o) : 


<«>rr #*, 


Tv/(-')»l,-l */<-«> 


u*(t) = - W; , 8 t P X + A x t ',f 

= - w;'b t p x + w;' b t x^*- a- x~’Pf 
= - w;' b t p x f w;' b t x,; t a _t xj p f 

= -w:'b t Px * w 3 'b t x* t a -t x* f 


“* » - Wj*' B T P X + Wj B A‘ T Pf 

Where A- T = -P(PA +W,)""'= (A-BW 3 'B T p) 

a(x II Ax;;)' T -x‘; ,, "A-'x 1 ‘; 0 

These are several forms for the control law. 
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1C) 


*(t)=(-Wj'B T p) x(t) 

+(wj'B T A' T p) f (+) 

where A=A“BW,B P 




limir 
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licqfjon to our 


roblem : 


U+)= C p x(t) +C v x(t) 


wh 


+ C d drt) + C . d (t) 

O w | 

cre ^J 0 » anc ^ ^<1, 

are constant chains. 
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Tal)lo 1 Optimal F/F and F/B Gains for Selected State 
Variable and Control Weightings. 


System Parameters: 

m = 100 lbm k = 0.3 Ibf/ft 

c = 0.000622 lbf-sec/ft (C = 0.17.) 
ft = 10 lbf/amp 


V <• 

ighls 


“la 

W ll> 

w 3 


92 6 


209 9 


258 

10 

581 

15 

1034 

21 ) 

1 (> 1 7 

25 

2320 

30 

1171 

35 

41 13 

10 

0325 

60 

1658 1 

80 

2591 1 

100 


F/B Chius 

0 (’ C r 

l> v l d0 'dl 

1.3845 1.3637 0.0294 - 0.0006 

3.1324 1.9863 0.0297 - 0.0001 

4.7659 2.4413 0.0298 - 0.0000 

6.3732 2.8210 0.0299 0.0000 

7.9701 3.1544 0.0299 0.0000 

9.5617 3.4552 0.0299 0.0000 

111 950 3 . 7354 0 . 0299 0 . 0000 

12.8153 3.9949 0.0299 0.0000 

I4 . 426' 1 4.2380 0.0299 0.0000 

I 

16.0321 4.4674 0.0299 0.0000 

24.0710 5.4729 0.0300 0.0001 

32.1259 6.3209 0.0300 0.0001 

10.1819 7.0680 00300 0.0001 

48.2297 7.7431 0.0300 0.0001 

56.2816 8.3640 0.(1300 0.0001 

64.3361 8.9420 0 . 03 ( H ) 0.0001 

96.5360 10.9526 0.0300 0.0001 

128.7372 12.6475 0.0300 0.0001 

160.9389 14.1407 0.0300 0.0001 


F/F Cains 

( '<12 C d3 

- 0.0070 - 0.0067 

- 0.0030 - 0.0019 

- 0.0020 - 0.0010 
- 0.0015 - 0.0007 

- 0.0012 - 0.0005 

- 0.0010 - 0.0004 
- 0.0008 - 0.0003 
- 0.0007 - 0.0002 

- 0.0006 - 0.0002 
- 0.0006 - 0.0002 
- 0.0004 - 0.0001 

- 0 . 0003 - 0.0001 

- 0.0002 - 0.0000 
4)0002 - 0.0000 

- 0.0002 - 0.0000 
4 ). 0001 4 ) . 0000 

- 4 ) 4)001 4 ) 4)000 

4 ) 4)001 - 0.0000 

-0 4)001 - 0 . 0000 


Ml 

-0 . 0049 
-0 . 0009 
-04)004 
- 0.0002 
-04)001 
- 0.0001 
-0.00(4 
- 0.0001 
- 0 . 0000 
- 0 . 0000 
-04)000 

-0 . oooo 
-0 . 0000 
- 0.0000 
-4)4)000 
4) . 00 JO 
-04)000 
-0 . 0000 
-0 0000 
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i nc 1 ude C bot h^Q^F/i'and propStioJal F/F;’C2 . 04, C6 and C8 
include LQR F/B alone. Weignung parameters used were 

_ A ^ A 1 / a /% i f\ r\ I 1 1 


“la = 258 - 


'lb 


j* a. v/ ** v r 

10, Wo = 1 (see Table 1). 


System Parameters 

Closed Loop Transfer Function 

s 2 X(s) . 

s 2 D( s) 

— i 

k (j^-) c ( lb f~— ) ■ C lbm) 

0 3 0.000622 100 

(C=0. 17.) 

0.3 0.000622 100 

0.0000622s +0.0001 
G 1 ( s ) - •? 

0 . 3 1056s “ +4 . 4675S+ 16 . 0o24 

0.0000622s + 0.0300 
G2(s) - 0 i 

0.31056s +4. 4675S+ 16. 0624 

0.45 0.000622 100 
0.45 0.000622 100 

0.0000622s + 0.0151 | 

b 0.31056 s 2 +4.4675s+16.0T74 ! 

0.0000622s + 0.0450 
S 0.31056s 2 +4.4675s+16.0774 

0.3 0.00622 100 

. 0.000622s + 0.0001 

° ^ 0.31056s 2 +4.4680s+16.0624 

0.3 0.00622 100 

0.000622s + 0.0300 j 

^ ~ 0 . 31056s 2 +4 . 4680S+16 . 0624 | 

0.45 0.00622 90 

0.45 0.00622 90 

.. . 0.000622s + 0.0151 

G((S) 7 , , 

0.27950s“+4.4680s+16.0( . 1 

0.000622s - 0.0450 j 

GS(s) - 7 — ~ " TI , 

0.27950s"+4.4680s+16.0i . 1 


ORIGINAL PAGE IS 
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Frequency (rad/sec) 


Figure 9 
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one 


I u s i 


ions: 


1. A n op + imal control has been determined 
for the nonkowogentous LQ R problem. 

2. An approximation to this optimal control 
has been -found which uses constant 
feedback and -feedforward gains. 

3. The optimal control has the following 
advantages: 

a. The gams can be easily determined. 

b. The control is very robust ^40*phase 
margin^ infinite positive gam margin 

4 dB negative gain margin). 

The control is applicable to a wide 
range of problems. 

The control offers substantial 
improvement in disturbance rejection 
over that afforded by LQ R feedback 
alone. 

e. The control can be easily implemented. 
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